Quantum computing with spatially delocalized qubits 
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We analyze the operation of quantum gates for neutral atoms with qubits that are delocalized 
in space, i.e., the computational basis states are defined by the presence of a neutral atom in the 
ground state of one out of two trapping potentials. The implementation of single qubit gates as well 
as a controlled phase gate between two qubits is discussed and explicit calculations are presented 
for rubidium atoms in optical microtraps. Furthermore, we show how multi-qubit highly entangled 
states can be created in this scheme. 
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Optical lattices and arrays of optical microtraps |^ 
are promising candidates for the implementation of quan- 
tum information processing with neutral atoms. Many 
of the requirements for quantum computation with opti- 
cal microtraps have been recently demonstrated in 2D 
arrays of ~ 80 traps with ^ 100 atoms per trap 
Some remarkable characteristics of optical microtraps 
are the possibility to scale, miniaturize and parallelize 
the required atom optics devices. In addition, they 
offer two fundamental advantages over optical lattices: 
(i) the possibility of individually addressing single traps 
due to the large separation of the microlenses foci, e.g., 

125/^m; and (ii) the independent displacement of rows 
and columns of microtraps and, eventually, of single mi- 
crotraps. Single atoms in dipole traps j^] and the Mott 
insulator transition with one atom per trap in optical lat- 
tices have been reported, and, therefore, the achieve- 
ment of ID and 2D arrays of optical microtraps contain- 
ing none or one atom per trap in a deterministic way can 
be foreseen for the near future. We will make use of all 
these features of optical microtraps to propose a novel 
implementation for quantum information processing. 

In our scheme, each qubit consists of two traps sep- 
arated by a distance 2a and one single atom. Per def- 
inition, the detection of the atom in the ground state 
of the left trap represents [0) and in the right trap |1), 
i.e., |0) = |0)^ and |1) = |0)^, where |0)^ are the 
vibrational ground states of the left and right trap, re- 
spectively. Throughout the paper we will call this im- 
plementation the spatially delocalized qubit (SDQ), since 
(0|r|0) — (l|r|l}| = 2a with r the position operator. To 
implement the SDQ we will assume that we are able to 
deterministically store none or one single atom per trap 
and cool it to the vibrational ground state in 3D. 

Single and two-qubit gate operations will be performed 
by adiabatically approaching two traps which will be 
modeled as follows: The initial separation of the traps 
is 2a„ia2;- The process of approaching them to the mini- 
mum separation 2amm takes a raising time tr. The tem- 
poral evolution of the distance a is described by the first 
half of a period of a cosine. The two wells remain at 
the minimum separation for an interaction time U and, 
finally, are adiabatically separated to the initial distance. 



To simplify the numerical analysis we will assume piece- 
wise harmonic trapping potentials as in ref. [|| and, even- 
tually, consider realistic Gaussian potentials as they are 
present in the experiment ^, ^ . 

Single-qubit operations, e.g., a Hadamard gate, are 
performed by adiabatically approaching the traps and 
allowing tunneling to take place. In order to illus- 
trate this operation, it is convenient to consider the 
two lowest energy eigenstates of the double well poten- 
tial. These two states are symmetric and antisymmet- 
ric, denoted by \S) and |^) respectively, with energies 
Es.Aia) — E{a) =F i?iSl(a) and H, being the splitting fre- 
quency. In terms of these states, our qubit basis reads: 
|0) = ^ i\S) + \A)), and |1) = ^ QS) - \A)). Let us 
assume that the atom is initially in the left trap, i.e, 
= 0)) = |0), then it is straightforward to check that 
its time evolution will be given by 



\m) 



cos( \Q)+ism(^t]\l) 



(1) 



Thus, the atomic wavefunction oscillates in a Rabi-type 
fashion between left and right traps at the flopping fre- 
quency iil(a). Obviously, for large trap separations 
states \S) and \A) become degenerate in energy, i.e., 
fl{a) = for a —>■ oo, and then the atom does not evolve 
in time (up to a trivial phase). Therefore, it is possible 
to realize single-qubit operations via tunneling by exper- 
imentally controlling the "Rabi frequency" through aniin, 
tr and ti. 

These single-qubit operations are illustrated in Fig. 1 
through a numerical integration of the ID Schrodinger 
equation in the parameter plane tr versus ti with |0) be- 
ing the initial state. Fig. 1 shows (a) the population 
of state |1), denoted by pi, and (b) the total ground 
state population of left and right traps, i.e., po + Pr, after 
the whole cycle of approaching and separating the traps. 
This oscillating population resembles the Rabi-flopping 
in the interaction of a laser with a two-level system and, 
for this reason, we have added in Fig. 1 dashed lines indi- 
cating the ^, TT, 27r laser pulse notation conventionally 
used in quantum optics. For small tr, non-adiabaticity 
results in the population of excited vibrational states, 
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FIG. 1: Single qubit operations for aamax = 5, aomm = 1-8 
and the atom initially in the left trap: (a) Population of the 
right trap; (b) The sum of the ground state population of left 
plus right traps. — ^Jhjnujj^ is the position spread of 

the ground state with m the mass of the neutral atom and l^Jx. 
the trapping frequency. 



which, as shown in Fig. 1 (b), yields po + Pi < 1- In 
what follows, wc will focus only in the adiabatic regime. 

For the two-qubit gate operations we assume that the 
two qubits are arranged either (i) in a ID configuration, 
i.e., the four traps form a line, or (ii) side-by-side in a 2D 
configuration, i.e., the traps form a square, c.f. Fig. 2(a). 
The traps involved are labeled Aq, A\ for the first qubit 
and Bo, -Bi for the second, and the respective ground 
states are denoted |0)yi, \\) a and |0)_b, A two- 

qubit gate will be realized via the coUisional interac- 
tion between bosonic neutral atoms 0. We will con- 
sider *^Rb atoms whose coUisional interaction can be de- 
scribed by a contact potential of the form Vifx^f'i) = 
47raf ?i^m~^(5'^(ri — r^), where at = 106 ao — 5.61 nm 
is the s-wave scattering length in the spin triplet. Be- 
cause the states |0) and |1) are localized in different po- 
sitions, it is enough to perform a suitable spatiotempo- 
ral variation of the potentials in order to pick up a col- 
lisional phase shift e.g. only if both atoms are in |1). 
This is sufficient to implement a two-qubit phase gate, 
which transforms product states \t)A\3)B-, hj G {0,1}, 
into exp{i6ii6jiTr)\i)A\j)B and, supplemented by arbi- 
trary single-qubit gates, forms a universal set of gates. 

For the in-line arrangement, the change of the poten- 
tial leading to a phase gate is shown in Fig. 2(b), where 



horizontal and vertical axes denote space and time, re- 
spectively. As the first step a tt pulse is applied on the 
second qubit, exchanging |0)b and \1)b- During this step 
only single-particle phases arise which can be included 
into the definition of the single-particle states. If the 
initial state was |0)yi|0)s then after the tt pulse, traps 
Ai and Bq would contain no atom. For initial states 
1)a|0)s or |0)yi|l)B, an atom would be either in Ai or 
in Bq and as seen before, we could approach and eventu- 
ally separate Ai and Bq such that a n27r pulse is applied 
with n integer. In this case initial and final state coin- 
cide, except for a single particle phase (ps which again 
can be included into the definition of |1)^ or |0)b. If we 
started from |1)a|1)b, then after the first pulse Ai and 
Bq would both be occupied, and during the n27r pulse 
the two atoms would collide. For an adiabatic evolution 
we can neglect not only the probability to populate ex- 
cited vibrational states, but also to find two atoms in the 
same trap, since, due to the coUisional interaction, these 
states are not degenerated with states where each atom 
occupies a different trap |5|. Thus for at initial and 
final state are the same except for a phase (j) and, in or- 
der to realize the desired phase gate operation, we need 
its coUisional part (j)c = 4> — 2<j)s to be an odd multiple 
of TT. The fidelity of this operation can be expressed as 
F — p - (cos(0c ~ "') + l)/2, which is plotted for the adia- 
batic regime in Fig. 2(c). Here p is the final probability to 
find the atom in the same state as it was before the n2'K 
pulse, neglecting the coUisional phase. To calculate the 
coUisional phase (j)c^ we have integrated the two-particle 
ID Schrodinger equation replacing t7(ri,r2) by an effec- 
tive ID interaction potential under the assumption that 
no transverse excitations occur Finally, to complete 
the phase gate operation another tt pulse is applied to 
the second qubit. 



In the case of 2D arrays of traps, as they are typically 
realized in the experiment Q| , the easiest operation is to 
move complete columns of microtraps. To realize the gate 
it is enough to be able to move selectively some columns, 
with the additional benefit that the operation is applied 
to many pairs of qubits in parallel which might allow for 
an easy implementation of error correcting codes. For 
the side- by-side arrangement, the initial and final tt pulse 
can be omitted and only the 27r pulse between traps Ai 
and Bi is needed. Although conceptually much more 
easier, the implementation in this arrangement demands 
the ability to move single traps instead of columns which 
makes it experimentally more involved. 



The SDQ scheme allows also for the creation of maxi- 
mally entangled states in a very straightforward way. Let 
us consider again four microtraps in the side-by-side ar- 
rangement. Fig. 3(a), with the two atoms located in the 
upper-left and lower- right traps, respectively, i.e., the ini- 
tial two particle state is |0)yii|l)B. During an approach of 
the traps. Fig. 3(b), the general two-particle state of the 
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FIG. 2: Two-qubit phase gate operation, (a) Arrangements 
of the qubits: (i) in-line, and (ii) side-by-side; (b) Contourplot 
of the spatiotemporal variation of the trapping potential for 
arrangement (i). The centers of the four traps are white, 
dark grey means high potential energy; (c) Fidelity F=p ■ 
(cos((^c — 1") + 1)/2. Parameters as in Fig. 1 with at = 106 ao, 
LLi^=ijy ^2 - and co^ = 1.1 ■ lO'^ s"\ 

system will be: 

= coo|0)a|0)s+coi|0)a|1)b + 
cio|1)a|0)b +cii|1)a|1)b + 

E c^io)^ii)j + E E (2) 

j=A,B i=0,lj=A,B 

Thus, the state of the system includes the four states of 
the computational basis, but also double qubit occupa- 
tion, defined as P dq=J2 j=A b '^j'^p ^^'^ double trap oc- 
cupation, pdt= J2i=o,i J2j=A,B ^ij^ij- order to create 
maximally entangled states, the approach of the traps 
must be both, adiabatic and symmetric, i.e, a{t) — b(t). 
Again, adiabaticity means that the population of excited 
vibrational states, as well as of double trap occupation 
states, can be neglected, c.f. Fig 3(c). A symmetric ap- 
proach results in poo{t) = Piiit) — /9dq(^)/2 during the 
whole process. In particular, these populations oscillate 
at the same frequency and, at the oscillation nodes, the 
state of the system is a combination only of |0)^|1)_b 
and |1)a|0)b. Therefore, choosing appropriate parame- 
ter values, such as those of Fig. 3, it is possible to obtain 
a maximally entangled state at the end of the process. 
Clearly, it is straightforward to generalize this approach 
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FIG. 3; Creation of a maximally entangled two-qubit state 
with two *^Rb atoms in a single step: (a) The initial state 
is |0)Ajl)s, with 2a and 26 being the separation of the traps 
in X and y direction, respectively; (b) We adiabatically and 
simultaneously, i.e., a(f) = 6(t), approach the four traps 
towards the center of the square; (c) Above: Population 
of states |0)a|1)s and |1)a|0)b; Below: Population of the 
state |0)yi|0)s (equal to the population of |1)a|1)b), dou- 
ble qubit, and double trap occupations. The parameters are: 
damage = 5, aami-n = 1.9, LUxU — 80, ujxti = 58, at = 106ao, 
ujx=Ljy ^2 - 10^ and = 1.1 ■ 10*^ 



to more than two qubits and, consequently, to prepare 
multiparticle entangled states in one single step. 

Although we have assumed harmonic trapping po- 
tentials so far, the experimental situation is de- 
scribed by Gaussian potentials of the form V{x) = 
— Vo exp(— (1/2) TOw^a;^/Vb). An analysis of the energy 
eigenvalues and eigenstates of the superposition of two 
such potentials as a function of the trap separation shows 
that the cosine function previously used to adiabatically 
approach the traps leads to values of tr larger by more 
than two orders of magnitude compared to harmonic 
traps. We have therefore applied the techniques from |^ 
to optimize the temporal variation of the trap separation 
while suppressing the population of excited vibrational 
states. For Vq = 200 huj^, Fig. 4(a) shows the result 
of this optimization for single-qubit operations, where 
the optimization is done with respect to the symmetric 
ground and first excited states. We notice that for the 
minimal distance 2amin the two traps are no longer sep- 
arated by a tunneling barrier, but they form a single fiat 
trap. Fig. 4(b) shows the population of the right trap, pi. 
The error rate due to the excitation of other vibrational 
states can be made smaller than 1% for w^tr > 1100 
which is a reduction by one order of magnitude compared 
to the non-optimized function a{t). 

Typical trapping frequencies for state-of-the-art 2D op- 
tical microtraps arrays of ®^Rb atoms are 10^-10^ s~^ in 
the transverse directions and IC'-IO^ s~^ along the laser 
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FIG. 4: Single-qubit gates for Gaussian trapping potentials 
for the following parameters: ujx ~ 6-10^ s~^, = 35.4 nm, 
Vo = 200hujj: = 0.9mK • ke, afflmai = 70, aa^in = 14.35: (a) 
The optimized variation of the distance for approaching the 
traps, insets show the form of the potential; (b) Population 
pi of the right trap showing Rabi-type oscillations. 



beam direction q| . In order not to populate excited vi- 
brational states and neglect double trap occupation, the 
trap displacement has to be adiabatic with respect to 
the lowest relevant trapping frequency. This adiabaticity 
condition yields realization times for single and two-qubit 
operations on the order of 1 ms and 10 ms, respectively. 
These values should be compared with the typical life- 
time of the atoms in the microtraps of 1 s, and the rate 
of spontaneous scattering of photons from the trapping 
laser of 1-10 s~^. The rate of decoherence for qubits en- 
coded in external states is further reduced by a factor of 
10-100 due the strong localization and almost harmonic 
trapping potential [nil . In addition, sideband cooling can 



be used to cool the atoms to the ground state of each trap 
in all dimensions. 

As discussed here, the SDQ configuration with neutral 
atoms in optical microtraps presents important practi- 
cal advantages with respect to the use of either inter- 
nal states [0, H or vibrational states ||, |l^ to encode 
the qubit. (i) Spontaneous emission leads to decoher- 
ence only in a much reduced fashion. As long as the 
microtraps are moved adiabatically, atoms remain in the 
ground state of both the internal and external degrees 
of freedom for typical gate times, (ii) There is no mo- 
mentum transfer in single or two-qubit operations as, in 
general, it is the case when these operations are realized 
via laser pulses. This momentum transfer could heat the 
atoms and, eventually, take them out from the microtrap. 
(iii) There is no need of state dependent interaction to 
realize the two-qubit operations. In the SDQ configu- 
ration one has that |(0|r|0) — (Ijr'll)! 7^ and, therefore, 
we can make use of the fact that all interactions are space 
dependent to realize the gate, (iv) The measurement of 
the state of the system is straightforward. Optical micro- 
traps can be separated to distances well beyond ~ 10/im 
which allows to detect the population of the trap by fo- 
cusing a laser field in one single trap and detecting the 
fluorescence signal Finally, (v) single and two-qubit 
gates are realized by the same kind of operation, i.e. by 
approaching the microtraps, which implies a strong sim- 
plification in the experimental set-up. 

Finally, it is important to note that most of the con- 
cepts developed in this Letter can be also applied to 
quantum dots with the qubit encoded in two tunnel-split 
ground states and the Coulomb interaction between elec- 
trons used to realize the qubit operations [|l^; and to 
Josephson-junctions based on the charge degree of free- 
dom with the Cooper-pairs tunneling coherently through 



the superconducting junction 1 13 . 
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